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Abstract 

In this paper we investigate the possible direct, non-gravitational interaction between holographic dark energy 
(HDE) and dark matter. Firstly, we start with two simple models with the interaction terms Q œ pg, and Q œ Page, and 
then we move on to the general form Q œ pi The cosmological constraints of the models are obtained from the 
joint analysis of the present Union2.1+BAO+CMB-+A data. We find that the data slightly favor an energy flow from 
dark matter to dark energy, although the original HDE model still lies in the 95.4% confidence level (CL) region. For 
all models we find c < 1 at the 95.4% CL. We show that compared with the cosmic expansion, the effect of interaction 
on the evolution of pam and pae is smaller, and the relative increment (decrement) amount of the energy in the dark 
matter component is constrained to be less than 9% (15%) at the 95.4% CL. By introducing the interaction, we find 
that even when c < 1 the big rip still can be avoided due to the existence of a de Sitter solution at z  —1. We show 
that this solution can not be accomplished in the two simple models, while for the general model such a solution can 


be achieved with a large £, and the big rip may be avoided at the 95.4% CL. 
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I. INTRODUCTION 


Cosmological observations such as the type Ia supernovae (SNIa) [1], the cosmic microwave background 
(CMB) [2] and the large scale structure (LSS) [3] all indicate that the universe is undergoing an accelerating 
expansion. This implies the existence of a mysterious component, named dark energy, which has negative 
pressure and takes the largest proportion of the total density in the present universe. In the last decade, lots 
of efforts have been made to understand dark energy [4], yet we still know little about its nature. 

In this paper we discuss the possible direct, non-gravitational interaction between the dark sectors in the 
framework of the holographic dark energy model, which is a quantum gravity approach to the dark energy 
problem [5-7]. In this model, the vacuum energy is viewed as dark energy, and is related to the event horizon 
of the universe when we require that the zero-point energy of the system should not exceed the mass of a 


black hole with the same size. In this way, we have the holographic dark energy (HDE) density [8] 
pae = 3C MR, (1) 


where c is a dimensionless model parameter, which can only be determined by observations, Mp; is the 
reduced Planck mass, and R; is the future event horizon of the universe, defined as 

Raf Caaf m Q) 
The HDE model has been proved to be a competitive and promising dark energy candidate. It can theoret- 
ically explain the coincidence problem [8], and is proven to be perturbational stable [9]. It is also favored 
by the observational data [10]. For more studies on the HDE model, see, e.g., [11-14]. 

The HDE model with some interaction between dark energy and dark matter (hereafter IHDE model) 
was firstly studied by Wang et al. in [15]. If dark energy interacts with cold dark matter,the continuity 
equations for them are pde + 3H(oqe + Pde) = —Q and Pam + 3Hpam = Q, where Q phenomenologically 
describes the interaction. The interaction between the dark sectors in the HDE model has been extensively 
studied in, e.g., [16]. It was found that the introduction of interaction may not only alleviate the cosmic 
coincidence problem, but also help to avoid the future big-rip singularity. 


There are various choices for the forms of Q. The most common choice is 


Q - FHp, (3) 


where T is a dimensionless constant, and p is taken to be the density of dark energy, dark matter, or the sum 
of them. These models are mathematically simple, so they are useful for phenomenology. However, it is 


difficult to see how they can emerge from a physical description of dark sector interaction. It is expected 


that the interaction is determined by the local properties of the dark sectors, 1.e., Pam and pae, but it is hard to 
understand why the interaction term must be proportional to the Hubble expansion rate H. A more natural 


form of the interaction was proposed by Ma et al. [17], where the interaction term takes the form 


Q « p$, pf. (4) 


In this description the interaction term only depends on the local energy densities of the dark sectors, and 
is thus more physically plausible. When a = 1 and f = O0, the interaction term Q œ pj; is the exact form 
we expected in the case of the dark matter decay. Similarly, the case of a = 0 and B = 1 corresponds to the 
dark energy decay. Moreover, in more complicated forms of interaction, like scattering, one may expect the 
existence of both pgm, and pae. Thus, Eq. (4) seems a more natural and physically plausible form to describe 
the interaction. 

In this paper we revisit the interaction between the HDE and dark matter by considering the interaction 
with the form Q « gu. To perform an overall analysis, we will take into consideration as many factors as 
possible. Different from many works in which the contribution of baryons is ignored, we will consider both 
baryon and dark matter components in the interacting models. Moreover, we will also discuss the models 
in a non-flat universe. As argued in Ref. [18], the studies of dark energy, and in particular, of observational 
data, should include Qyo as a free parameter to be fitted alongside the w(z) parameters. Another reason why 
we take the spatial curvature into consideration is the possible correlation between the curvature and the 
interaction. As pointed out in Ref. [19], compared with the flat universe, a much stronger interaction can be 
allowed by the data in a non-flat universe. 

There are many interesting issues worth investigating in the IHDE models. For example, the properties 
of the interaction term, the direction of the energy flow, the evolution of the dark matter and dark energy 
densities, the fate of the universe, and so on. In the IHDE model with Q « pe. ph. some of these issues have 
been discussed in [17]. In this paper, we will perform a more comprehensive exploration of these issues. 

This paper is organized as follows. In Sec. II, we derive the basic equations for the IHDE models. In 
Sec. IIL, we introduce the methodology and data used in this work. In Sec. IV, we discuss the cosmological 
interpretations of the three IHDE models, including two simple models with Q œ pgm, Q « Pde and one 
general model with a and f treated as free parameters. We firstly study the cosmological constraints on these 
models, and then discuss the fate of the universe in the models. At last, we give some concluding remarks 
in Sec. V. In this work, we assume today's scale factor ag = 1, so the redshift z satisfies z = 1/a — 1; the 
subscript “0” always indicates the present value of the corresponding quantity, and the unit with c = A = 1 


is used. 


II. INTERACTING HOLOGRAPHIC DARK ENERGY MODEL IN A NON-FLAT UNIVERSE 


In this section, we give the basic equations for the IHDE model in a non-flat universe. 


A. Friedmann equations in a non-flat universe 


In a spatially non-flat Friedmann-Robertson-Walker universe, the Friedmann equation can be written as 


3M;H* = Pdm + pb + Pr + Pk + pde: (5) 


where px = -3M2,5 is the effective energy density of the curvature component. For convenience, we define 


the fractional energy densities of the various components, i.e., 


-k 
= P Qde = LA Quam = ee > Q, Tea Q, = E (6) 


E- Pe Pe Pe Pe Pe 


Qg 
where pe = 3M? ,H? is the critical density of the universe. The subscripts, k, de, dm, b and r, represent 


curvature, dark energy, dark matter, baryon and radiation, respectively. By definition, we have 
Qde F Quam F Q, F Q, =F Ox =1. (7) 


With the existence of interaction between the dark sectors, the energy conservation equations for the 


components in the universe take the forms 


Pam + 3Hpam = Q, (8) 
Pde + 3H (Pde + Pade) = —Q. (9) 
Dp + 3Hpy = 0, (10) 

pr + 4Hp, = 0, (11) 

pe + 2Hpy = 0. (12) 


As mentioned above, Q denotes the phenomenological interaction term. Combining Eqs. (8)-(12) together, 


we can obtain the form of pae, 


2H 1 1 
Pde = — 3 gabe Pe apr + ape (13) 


Substituting Pae into Eq. (9), we have 
QË + Oe 
H Que 


+ 3H pae + H(px — pr) — yf +3H Pc --Q. 
H 
Dividing the above equation by pe, we get a derivative equation of H and Qye, 


H . 
2(Que > Dg T Que T H(3904, -3+ [971 x Q,) = —HXj, 


where we have defined the effective dimensionless quantity for interaction, Qz, with the form 


Q 


Q; = ——. 
H(z)pc 


B. HDE in anon-flat universe 


From the energy density of the HDE, Eq. (1), we have 


C 


L= = 
H V Oe 


Following Ref. [13], in a non-flat universe, the IR cut-off length scale L takes the form 
L — ar(t), 


and r(t) satisfies 


f” dr {~ dt 
0 V1 — kr? 7 t a(t) 
By carrying out the integration, we have 


] , TE ] + da 
r(t) = ment [ 2! = TE sin ( Vk T 22) 


Equation (18) leads to another equation about r(t), namely, 


C 


L 
(n = — = ——_.. 
- a VQaeHa 


Combining Eqs. (20) and (21) yields 


+œ dt k 
Vk f — = arcsin alae 
t a VOgeaH 


Taking derivative of Eq. (22) with respect to t, one can get 


Ò 1 [oum 
de +H + z = Qa H^ — k 
2Qde H c? a? 


(14) 


(15) 


(16) 


(7) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


C. Evolution equations of E(z) and Q4e(z) in an IHDE scenario 


Combining Eq. (15) with Eq. (23), we eventually obtain the following two equations governing the 


dynamical evolution of the IHDE model in a non-flat universe, 


1 dEG) — Og (Q-Q,-3*0; 1, [Ou |, (24) 
E@ dz 1+z 208, f ET 


Adde _ eo We i ù%-Q,+Q; 
Se ee a PAg —— — —— |. 


2. 20- Qn) n 


dz 1 +z 
where E(z) = H(z)/ Ho is the dimensionless Hubble expansion rate. Equations (24) and (25) can be solved 
numerically and will be used in the data analysis procedure. Notice that we have 


Oyo(1 +2)? Q,o(1 + z)* Opd + 2° 


eg egg. SUN pay 


à (26) 


and the fractional density of dark matter is given by Qam(z) = 1 - Q(z) — Qae(z) - Q(z) - Qp(z). The values 


of O5 and Oo, for simplicity, are determined from the 7-yr WMAP observations [23], 


Quo = 0.02253h 7, (27) 


Qro = Qyo(1 + 0.2271 Ners), Qo = 2469 x 10?n 2, Neff = 3.04, (28) 


where y represents photons, and Nepp is the effective number of neutrino species. 


D. Models 


In this paper we consider the interaction term with the form Q o ph. For concreteness, we express 


the interaction term as 


a 
p dm i 


Q =T Ho -pT (29) 


cO 
which can also be expressed as Q — TH op co F(z) P Qam) Qef, and the parameter I’ is dimensionless. 
We will investigate three IHDE models in this paper. Firstly, we discuss two simple models with fixed 


a and f. The first model (hereafter IHDE1) is the case of a = 1 and £ = 0, and thus 
Q -FHopas, Qi = TQam(Z)/E@). (30) 


When [ > 0, the energy transfer corresponds to the decay of dark energy into dark matter; vice versa. 
This model is similar to the model proposed in [20], where the same form is introduced in the interaction 


between dark matter and quintessence. 


The second model (hereafter IHDE2) considered in this paper is the case of a = 0 and B = 1. Corre- 


spondingly, we have 
Q=THopae, Qi = TQae(Z)/E@). (31) 


Finally, we also consider a more general model (here after IHDE3) where a and f are treated as free 
parameters. The formula of Q has been given in Eq. (29), while the dimensionless quantity for interaction 


defined previously is 


Or = TEGO? Qam) Oy. (32) 


III. METHODOLOGY 


For the IHDE models in a non-flat universe, there are six free parameters: c, O50, Quo, T, œ and f£. We 
will constrain them by using the latest observational data. As a comparison, the ACDM model and the HDE 
model with spatial curvature but without interaction (namely, €2,9 + 0 but Q = 0) will also be investigated. 

In this work, we adopt the y statistic to estimate the model parameters. For a physical quantity £ 
with experimentally measured value ops, standard deviation cz and theoretically predicted value £j, the x 


function takes the form 


2 (Eobs = En)? 


= M 33 
XE = (33) 
é 
The total y? is the sum of all AGS: Le. 
XY». (34) 
é 


One can determine the best-fit model parameters by minimizing the total y*. Moreover, by calculating 
Ay a7 = a us one can determine the 68.3% and the 95.4% CL ranges of a specific model. 

In this work, we determine the best-fit values and the 68.3% and 95.4% CL ranges of the model pa- 
rameters by using the Markov Chain Monte Carlo (MCMC) technique. We modify the publicly available 
CosmoMC package [21] and generate 109—107 samples for each set of results presented in this paper. 

For data, we use the Union2.1 SNIa sample [22], the CMB anisotropy data from the 7-yr WMAP ob- 
servations [23], the BAO results from the SDSS DR7 [24], 6dFGS [25] and WiggleZ Dark Energy Survey 
[26], and the Hubble constant measurement from the WFC3 on the HST [27]. In the following, we briefly 


describe how these data are included into the y? analysis. 


A. The SNIa data 


First we start with the SNIa observations. We use the latest Union2.1 sample including 580 SNIa that 
are given in terms of the distance modulus Hops(zi) [22]. The theoretical distance modulus is defined as 
Uthi) = 510810 Dr(Gi) + uo. (35) 


where uo = 42.38 — 5 log), A with h the Hubble constant Ho in units of 100 km/s/Mpc, and the Hubble-free 


luminosity distance Dz = Hod; is 


IQ 
D(z) = "sim vI9o| JE EG 5) (36) 
where 
sin(x), if Quo < 0, 
sinn (x) = 4 x, if O49 = 0, 


sinh(x), if Oo > 0. 
The x? function for the SNIa data is 


Y [Mobs (i) — pap (37) 


c? 
i 


X$N = 


where Hopbs(zi) and g; are the observed value and the corresponding 68.3% error of distance modulus for each 
supernova, respectively. For convenience, people often analytically marginalize the nuisance parameter uo 
(i.e., the Hubble constant Ho) when calculating XS n [28]. 

It should be stressed that Eq. (37) only considers the statistical errors from SNIa, and ignores the sys- 
tematic errors from SNIa. To include the effect of systematic errors into our analysis, we will follow the 
prescription for using the Union2.1 compilation provided in [29]. The key of this prescription is a 580 x 580 
covariance matrix, Cs y, which captures the systematic errors from SNIa (This covariance matrix with sys- 


tematics can be downloaded from [29]). Utilizing Cs y. we can calculate the following quantities 


A= (u — IP Cs Yit = in» (38) 
580 
B=) (Cyn i - um), (39) 
i=l 
580 
C=) Gin (40) 
i,j-l 


and the y? function for the SNIa data is [29] 


2 


2 
mc 4l 
XSN C (41) 


Different from Eq. (37), this formula includes the effect of systematic errors from SNIa. 


B. The CMB data 


Here we use the “WMAP distance priors” given by the 7-yr WMAP observations [23]. The distance 
priors include the "acoustic scale” /4, the "shift parameter" R, and the redshift of the decoupling epoch of 
photons z+. The acoustic scale /4, which represents the CMB multipole corresponding to the location of the 


acoustic peak, is defined as [23] 


lim (ez. (42) 


Here Da(z) is the proper angular diameter distance, given by 
Dalz) = di (Of +2)", (43) 


and r,(z) is the comoving sound horizon size, given by 


1 1/(1+z) da 
rz) = —= | es (44) 
V3 Jo a2H(a) 41 + 8059/4,9)a 
where (2,9 and Qo are the present baryon and photon density parameters, respectively. As mentioned above, 
we adopt the best-fit values, Qy9 = 0.02253h~? and Qo = 2.460 x 107547? (for Temp = 2.725 K), given by 


the 7-yr WMAP observations [23]. The fitting function of z. was proposed by Hu and Sugiyama [30]: 
z, = 1048[1 + 0.00124(Q59/7) 9511 + ¢1(Qnoh?)] , (45) 


where 


|. 0.0783(Quoh?) 0258 E 0.560 
-1«395(0,)076 ^. 7 TXZZLKQO)U | 


Here the subscript m denotes the matter component, i.e., Qm = Qam + Qp. In addition, the shift parameter R 


81 (46) 


is defined as [31] 
R(z.) = JOgoH?(1 + z.)D(.) . (47) 


This parameter has been widely used to constrain various cosmological models [32]. 


As shown in [23], the y? function of the CMB data is 
Xene = 8” — (Comp dis - XD. (48) 
where x; = (l4, R, zx) is a vector, and (Coup jis the inverse covariance matrix. The 7-yr WMAP observa- 
tions [23] have provided the maximum likelihood values: /4(z.) = 302.09, R(z.) = 1.725, and z. = 1091.3. 
The inverse covariance matrix was also given in [23], 
2.305 29.608 -1.333 


(Cclg) =| 29.608 682527 -113.180 |. (49) 
~1.333 -113.180 3.414 


C. The BAO data 


In this paper we use the BAO data from the SDSS DR7 [24], the 6dFGS [25], and the WiggleZ Dark 
Energy Survey [26]. In the following, we will describe the BAO distance measurements of these projects, 
and introduce how to add them into the X? statistics. 


One effective distance measure is Dy(z), which can be obtained from the spherical average [33] 
1/3 
Dy) =|(1 + 9° Di) — (50) 
A A H , 

where D(z) is the proper angular diameter distance. 

For the 6dFGS and SDSS DR7 data, we use the quantity d, = r;(zq)/Dy(z). The expression of r, is given 
in Eq.(44), and zg denotes the redshift of the drag epoch, whose fitting formula is proposed by Eisenstein 
and Hu [34], 


129 1(Q,,9//2)025! 


= | 1 bi (9:982)? | , " 
S 706550 gO | + bi (Quoh") | (51) 
where 
by = 0.313(Qn0h?Y ^? [1 + 0.607 Qnoh?)™], (52) 
bz = 0.238(Q,,952)9223. (53) 


For the SDSS DR7 data [24], we write Y for the BAO data as 
2 -1 
Xpao,s pss = AP(Cpao,s pss)ijAPj. (54) 
where 
Ap; = py" — pi, pi = dy" = 0.1905, pj** = dy45 = 0.1097, (55) 


and the inverse covariance matrix takes the form 


x 30124 —17227 
(Cp40,s pss) = (56) 
—-17227 86977 
The 6dFGS survey [25] gives a measurement of quu. = 0.336 + 0.015, so we have 
do.106 — 0.336 V 
XBAO6dFGS = M | (57) 


The WiggleZ Dark Energy Survey [26] gives three measurements of the A parameter in different red- 
shifts. The A parameter is defined by [33] 


A(z) = 


100Dy(z) YQnoh? (58) 
: 


10 


and the x? of the WiggleZ Dark Energy Survey takes the form 


2 _ -1 
X BAO,Wigglez = APi(CBA0,Wigglez)ijÂP j» (59) 
where 


Api = pi — pi, pi = Agta = 0.474, pi = AFE = 0.442, p^ = AUS = 0.424, (60) 


I 
and the inverse covariance matrix takes the form 


1040.3 -807.5 336.8 
(Crs wee) = | 8055 - 37208. »i5515]. (61) 
3368 -1551.9 2914.9 


The final BAO y is a combination of the SDSS DR7, the 6dFGS and the WiggleZ BAO ys Le: 


2 _2 2 2 
XBAO 7 XBAO,SDSS + XBAO6dFGS * XBAO,WiggleZ* (62) 
D. The Hubble constant data 


The precise measurements of Ho will be helpful to break the degeneracy between it and dark energy 
parameters [35]. When combined with the CMB measurement, it can lead to precise measure of the dark 
energy equation of state (EOS), w [36]. Recently, using the WFC3 on the HST, Riess et al. obtained an 


accurate determination of the Hubble constant [27], 
Ho = 73.8 + 2.4 km/s/Mpc, (63) 


corresponding to a 3.396 uncertainty. So the y? of the Hubble constant measurement is 


2 
" ( > e) 


é 4 
Xn 0.024 (6) 


E. The total y? 


Since the SNIa, CMB, BAO and Hp are effectively independent measurements, we can combine them 


by simply adding together the y? functions, i.e., 


2 2 2 2 2 
Xtotal ^ XSN * XcMB * XBAO + Xp (65) 
IV. RESULTS 


In this section, firstly, we show the cosmological constraints on the three IHDE models, and then we 


discuss the cosmological implications of these models. 
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TABLE I: Fitting results of the models 


Model Oo c c (95.496 range) Tr Quo Xx 

ACDM 0235570 - - O (fixed) ^ 0.001*0003 ^ 550.354 
HDE (Q - 0) 0.2233*0012  0.71'009 | 0.56« cx 0.94 — O(fixed) — 0.006'00077 549.461 
Q = FHopan 0.238+9013  0.69+010  0.54« c «0.91. —0.056:0031 0.015+0.011 548.352 
Q =THopae 0.23770-013 0.664011  0,50«c«0.90 -0.07340079 0.016002 548.390 
Q =H p [pt 0.239+9013 — 0,63*0.8  047<c<0.92  -0.014*0914  0.017+%011 548.298 


A. Cosmological constraints 


In this subsection, we discuss the cosmological constraints on the three IHDE models. The parameter 
space of these IHDE models are explored by using the Union2.1+BAO+CMB-+AH data. We summarize the 
fitting results in Table I. For comparison, the fitting results of the ACDM and HDE models by using the 
same set of data are also shown. 


The values of x? .. of the models are listed in column 7. Clearly, compared with the ACDM and HDE 


min 


2 


models, the IHDE models lead to a evident reduction of y; . 


However, taking the number of parameters 
into account, ACDM and HDE models still provide nice fits to the data. In addition, compared with the two 
simple IHDE models, the complicated model with Q o p. does not lead to a significant reduction of 
AN sd implying that in the context of the current observations, such a complicated model is not necessary. 

In column 2-6, the best-fit values and the 68.396 CL (A? = 1) uncertainties of the parameters Quo, c, 
T and Qy are listed. For the HDE and three IHDE models, we found c < 1 at the 95.4% CL (Ay? = 4). 
This result is consistent with the results in the previous works [10, 16]. For the HDE model, c < 1 means 
that the universe will end up with a the big rip, while for the IHDE models, due to the interaction between 
the dark sectors, the big rip may be avoided. We will discuss this issue in the following content. 

For these three IHDE models, we find T < 0 at the 68.3% CL, namely, the energy flow from dark matter 
to dark energy is slightly favored by the data. 

In the left panel of Fig. 1, we plot the contours of the 68.396 and 95.496 CL for the three IHDE models in 
the I—c plane. We can see that the upper bounds on the parameter c are similar, while the lower constraints 
are slightly different from each other. In the IHDE3 model, due to the complexity of the model, the allowed 
parameter space is much larger than those of the IHDE1 and IHDE2 models. 

Interestingly, it seems that, compared with the HDE model, the IHDE models slightly favor smaller 


values of c. This can be seen in the column 4 of Table I. To see it more clearly, we also plot the probability 
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FIG. 1: Marginalized probability contours at the 68.3% (Ay? = 2.3) and 95.4% (Ay? = 6.18) CLs in the I—c (for 
the three IHDE models) and Q4mo—c (for the HDE model) planes. The contours for the IHDE1 (cygan, labeled as I), 
IHDE2 (orange, labeled as II), and IHDE3 (olive, labeled as IIT) models are plotted in the left panel, and the contours 
for the HDE model (magnetic) are plotted in the right panel. The green dashed line denotes c = 1. Clearly, compared 


with the HDE model, the fitting results of c in the three IHDE models are smaller. 


contours in the Qgmo—c plane for the HDE model in the right panel of Fig. 1. The 95.4% CL of the HDE 
model intersects with the line c = 1, while the contours of the IHDE models are all below the line. 

In the left panel of Fig. 2, we plot the contours of the 68.3% and 95.4% CL for the three IHDE models 
in the I-O;9 plane. The parameter space of the IHDE3 model is much larger than those of the IHDEI and 
IHDE2 models, and we can see clear degeneracy between Qgo and I’, consistent with the result of [19]. The 
degeneracy amplifies the range of Qz9 compared with the HDE model without interaction, as is seen in the 
right panel of Fig. 2, where the contours in the 07,,5-O;9 plane for the three IHDE and the HDE models are 
all plotted. 

We are also interested in the constraints on the parameters a and f, i.e., the interaction forms allowed by 


the data. In our analysis, we find the 95.4% CL constraints are 
ge. B= 10.981387. (66) 


So the allowed regions of œ and f are very wide. This phenomenon is understandable: the intensity of the 
interaction term Q is almost decided by T. If T is small enough, any large values of œ and £ are tolerable. ! 


The 68.3% and 95.4% CL contours in the a—6 plane are plotted in Fig. 3, and we can see strong degeneracy 


! In our analysis we generate about 10’ samples for the IHDE3 model. It is expected that the area of the contour could be larger 
if one performs a more extensive analysis. 
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FIG. 2: Marginalized probability contours at the 68.3% and 95.4% CLs in the I-O;o (left panel) and Qumo—Qyo (right 
panel) planes for the IHDEI (cygan, labeled as I), IHDE2 (orange, labeled as II), IHDE3 (olive, labeled as III), and 
the HDE (magnetic) models. The Qy,9—Qyo contours of the three IHDE models are similar to each other, so for the 


IHDE2 and IHDE3 models, we only plot the 95.4% CL contours. 
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FIG. 3: Marginalized probability contours at the 68.3% and 95.4% CLs in the a— planes for the IHDE3 model. 


between a and £. 


B. Constraints on the interaction term 


In this subsection we discuss the constraints on the interaction term in the IHDE models. The fitting 


results of I have been listed in the 5-th column of Table I. We have seen that for all the models we have 


0.28 


Q « PH Pam pac’ 


FIG. 4: Reconstructed evolutions of 3H pam, 3H(1 + weff)Pae and Q (all divided by pco) at the 95.4% CL, for the 
IHDEI (left panel), IHDE2 (middle panel) and IHDE3 (right panel) models, respectively. 


T < Oat the 68.3% CL. Here we list the 95.4% CL constraints for the three IHDE models, 


-0.154 < T € 0.053 for IHDEI, (67) 
-0.229 < T < 0.059 for IHDE2, (68) 
-0.509 < T < 0.072 for IHDE3. (69) 


For the I > 0 case, i.e., the direction of the energy flow is from dark energy to dark matter, we find a tight 
constraint |I| < 0.07 at the 95.4% CL for all the three IHDE models. For the F < 0 case, the allowed value 
of T is much larger. At the 95.4% CL, for the IHDE1 and IHDE2 models we have |I] S 0.2, while for the 
IHDE3 model we have |I] < 0.5. 

To see the result more vividly, in Fig. 4 we reconstruct the evolution of 3H pgm, 3H(1 + Wae)Pde, and Q at 
the 95.4% CL for the three IHDE models. As shown in Eqs. (8) and (9), the term 3H pgm [or, 3H. +Wae) Pde] 
characterizes the influence of the cosmic expansion on the energy density of dark matter (or, dark energy), 
while Q describes the effect of the interaction. As shown in Fig. 4, for all the models Q is much smaller 
compared with 3Hp;,, and 3H(1 + Wae)Pae, implying that the influence of interaction on pam and bde is 
much weaker than the cosmic expansion. It is also interesting to find that the evolution of the interaction 
terms have different behaviors in the three IHDE models, due to the different expressions of Q. In the 
IHDEI model with Q œ pgm, the absolute value of Q may significantly decrease (along with pam) during 
the evolution, while in the IHDE2 model the absolute value of Q almost maintains no changes during the 
evolution since Q œ pge. 


We are interested in the influence of the interaction on the dark matter density, and thus we investigate 
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the quantity Q/(3Hpam), which describes the ratio between the changes of pg, caused by the interaction 


and by the cosmic expansion. Here we list the constraints on this quantity at z = 0 at the 95.4% CL, 


— 0.0514 < Q < 0.0177 for IHDE1, (70) 
Pdm 

-0.2097 < —— < 0.0625 for IHDE2, (11) 
Pdm 

—0.4779 < -— €0.0105 for IHDE3. (72) 
dm 


We see that in the IHDE1 model, the present effect of Q is much smaller than the other two IHDE models. 
The large value (47.79%) of the lower bound of this ratio in the IHDE3 model indicates that in this model 
the influence of Q on 6g, can be comparable with that of the cosmic expansion. Figure 4 shows that in the 
IHDEI model the present value of Q is much tightly constrained compared with the other two models, and 
in the IHDE3 model Q could take large negative values when z — 0. 

Also, since we are interested in how much dark matter is produced or annihilated due to the interaction, 


it is helpful to define another quantity, 


3 
Pdma 


m 
A= z=0 ` Pdma ki 


- (73) 
Pama b 
Notice that pama? is the total energy of dark matter in a unit comoving volume. In the absence of interaction, 


Pama? is conserved, and A = 0. Here Pink is the "initial" value of the total energy of dark matter in a 


ini 
comoving volumn. It can be calculated at the early epoch, e.g., z > 5000. Thus, A describes the relative 
change amount of energy in the dark matter component caused by the interaction. 


The results for the three IHDE models are as follows: 


— 15.22% < Ax 8.8396 for IHDEI, (74) 
-13.3% < Ax 5.5796. for IHDE2, (75) 
~14.49% < Ax 6.5796 for IHDE3. (76) 


We find that the results of the three models are not much different from each other. Thus, although the 
evolution of the interaction term is evidently different, the constraint on the overall effect of the interaction 
does not change a lot. Roughly, at the 95.4% CL, the increment amount of the energy of dark matter is 


constrained to be less than 9%, while the decrement amount is constrained to be less than 15%. 


C. The fate of the universe 


In this subsection we discuss the fate of the universe in the IHDE models. The most important reason for 


the introduction of the interaction between the dark sectors in the HDE model is to avoid the future big-rip 
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singularity. For the IHDE model with Q « prar. this issue has been briefly discussed in [17], where it 
is shown that back reaction effect from quantum corrections cannot prohibit the occurrence of the big rip. 
Here, we will give a more detailed investigation of the fate of the universe in the IHDE models considered 
above. 

Firstly, we will derive some useful formulas as a preparation of our discussion. Then, we will analyze 
the dynamical equations to see whether it is possible to avoid the big rip when c « 1. Finally, based on the 
constraints of the data, we will numerically solve the equations and reconstruct the evolution of the dark 


sectors to z — - 1. 


1. The effective EOS of dark sectors 


We can define the effective pressure and effective EOS of dark energy, i.e., 


Q 
Pef fe = Pde + aH Wef f,de = Wde + 3Hpae’ (77) 


satisfying 
Pae + 3H (Pade + Peff) =0, pae +3H(1 + Wef f,de)Pde = 0. (78) 


From Eq. (13), using 4 =-(1 *2Hf£ and E =-(1+ zH =-(1+ 2255, we derive the equation 


of state for dark energy, 


20-2dlE 1 19, 1 


= —— — - -— + -—. 79 
Me "30% dz Qu 304 30 va 
Combined with Eq. (24), it follows that 
1 2 Qz Q; 
Ta aii = , 80 
Wde 3 E. tS 304. (80) 
which reduces to the familiar formula wg, = -i — i Ox when Q; = 0 and Q; = 0. Thus, the effective EOS 


takes the form 


Weff,de = Wde + Te (81) 
1 2 [Q4¢e 

= -l2-l4[|—-a9, 82 

3 3V2 t% (82) 


This result is very interesting: We¢fde does not explicitly contain the interaction term. However, we should 
keep in mind that the evolutions of Qge(z) and E(z) are determined by the differential equations (24) and 


(25) where the interaction Q (or, Q7) is involved. 
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It is much easier to derive the effective EOS of dark matter. From Eq. (9), it follows directly that 


i Q | Qj 
eines 3H pan i 3045, ES 
- -3TEQ 93g 1of (84) 


To investigate the fate of the universe, we are also interested in the effective EOS of the total components 


in the universe, i.e., 


Wef ftot = P3 Wef fii, (85) 
j 


where i represents de, dm, b, r and k. In the future epoch, only the dark energy and dark matter components 


are important. 


2. Analysis of the dynamical equations in the z —^ —1 region 


Now we make some general investigations to see whether the big rip can be avoided when c « 1 in the 
IHDE models. Equation (82) shows that the effective EOS of HDE is only determined by O;, and c (we 
neglect Qg which is less important in the future). In the case of c < 1, once O,, — 1 there will always 
be Werfde < —l and Wer ftor < —1, and the fully dominated dark energy will drive the universe to a big 
rip. Thus, to avoid the big rip, we must evade O4, — 1 in the future. Also, to avoid pge from infinitely 
increasing, we must require I > 0, i.e., the direction of energy flow is from dark energy to dark matter. 

Here we list the simplified equations of motion in the region z  —1. Neglecting the radiation, baryon 


and curvature components, we have 


1 dE(z) — Qde [Qde -3 + Oj " VO, (86) 
E(z dz 1+z 20, c JJ 
dO, EN 2041 = Qae) VQde l = OQ, (87) 
dz E. c 2. 20—05,/. 


Notice that unlike the cases of Q œ Hpam or Q œ Hpge, the two equations are coupled, and thus the situation 


is much more complicated. However, we find that these equations have a stable point 
Que =, Qj 230 -= c°), (88) 


which leads to 


Weff,de = E E =-l, Wef f,dm == 3n. = zk (89) 


E T Ei T ST F T T T T 


| c-0.9, P=0.1, a=0, B=10 
| 


FIG. 5: Evolutions of E, O,, Qam, Oy, Weff,de and Weff am along with the scale factor a. Here we take typical values 


of Qang = 0.24, Ojo = 0, A = 0.73, and others are denoted in the panel. 


and 


dE@) dO 
dz” dz 


=0. (90) 


This is exactly a de Sitter solution. See also [37] for a relevant study about the avoidance of big rip in the 
holographic dark energy model, where the extra dimension effect is considered and the final state of the 
universe is also a de Sitter solution. 

We find that this situation can really happen when £ is large, and the interaction is strong enough to pull 
down €, from 1. As an example, in Fig. 5 we demonstrate how the quantities evolve along with the scale 
factor a in the case of B = 10. A sophisticated discussion on the evolution behaviors of the quantities in the 
whole range of a would be tedious and unnecessary. Thus, instead, we just focus on the epoch of a > 10 
(the region to the right of the thick dashed line). This would be enough for us to understand how it happens. 
In this epoch, firstly, E(z) increases due to the dominated dark energy component with weffde < —1, while 
the strong interaction pulls Qg. down from 1. So along with a, Qae decreases and correspondingly Qam 
increases. As a result, We fae increases from values less than —1, while the effective EOS of dark matter, 
Weffdm = a. which is highly negative due to a small Qy,, also goes up. When both wae eff and Wam,ef f 
approach —1, the system reaches the stable point, and E(z) stops increasing. The final state of the universe 
is de Sitter-like. 


Note that an essential condition that such a procedure can happen is that the interaction is strong enough 


to pull O4, down from 1. Thus, a large 6 is expected. To see the dependence on the values of f, in Fig. 6 
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FIG. 6: Evolutions of Q;, (left panel), Que (left panel), weff,de (middle panel) and Wes fam (right panel) along with the 


scale factor a for three sets of œ and £. The same as Fig. 5, we take values of O,,9 = 0.24, O,9 = 0 and h = 0.73. 


we plot the evolution of Que, Qdm, Wef fide and Wef f,dm in large a regions for different values of œ and B. For 
the case of a = 0 and £ = 1 (the blue line), the interaction is not strong enough to avoid Que 1, so the big 
rip still happens. For the case of a = 0 and B = 3 (the green dotted line), the interaction is strong enough, 
and so we see that Qge is pulled down from 1, and the system evolves to a stable point when a > 10000. 
For the case of a = —1 and £ = 7 (the red dashed line), the interaction is so strong that the system quickly 
approaches the stable point at a « 1000. 

In the following, we will numerically solve the equations for the IHDE models in a large a region in 
the parameter space constrained by the data, and see whether the de Sitter solution can be achieved in these 


models. 


3. TheIHDE! and IHDE2 models 


In the top panels of Fig. 7 we show the reconstructed evolution of Quy. along with z for the IHDEI and 
IHDE2 models. Clearly, at the 95.4% CL we see Qa. — 1 when z > —1, and we have weff,de < —1 since 
c « 1. So these models will not help us to avoid the big rip. 


This result is understandable. The fate of the universe in the IHDE model with the interaction term 
Q = 3bH pae (91) 
was discussed in [9] and it was shown that the condition to prevent the big rip is 


but", (92) 
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FIG. 7: The evolutions of Qu. (top panels) and wer fe (bottom panels) along with the redshift z at the 68.3% and 
95.4% CLs, for the IHDEI (left panels) and IHDE2 (right panels) models. The 95.4% CL of HDE model is plotted in 


dashed lines for a comparison. 


In the IHDE1 and IHDE2 models, analogously, we define the "effective" coupling, 


Q 
e = , 93 
p 3Hpae 
and we have 
IO, 
b, , berfp = 94 
^f^ 3E@Q4ue H ^ SEO id 


for the two models, respectively. Note that once E(z) becomes increasing in the future, besp will be sup- 
pressed to zero, so the condition (92) cannot be satisfied. 

The reconstructed wef f, de(z) for the two models are shown in the lower panels of Fig. 7. As a comparison, 
the 95.4% CL evolution of the HDE model is also plotted (the black dashed line). We see that for these 
models Wer fde(Z — —1) < —1, and thus the big rip will happen. Also, the values of weff,de in the IHDE 
models are more negative than that of the HDE model, since in the two IHDE models c is constrained to 
be smaller values. Thus, instead of helping us to avoid the big rip, the situation is even worse in these two 


models. 
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FIG. 8: The evolutions of Qae (gray, left panel), Qam (green, left panel), werfae (gray, right panel) and wef ftor (orange, 
right panel) along with the scale factor a for the IHDE3 model. The 68.3% and 95.4% CLs are plotted. 


4. The IHDE3 model 


We have seen in Fig. 6 that a large f can drive the system to a stable point and avoid the big rip. As 
shown in Eq. (66) and Fig. 3, in our fitting results we found that large values of B are allowed, so we expect 
that the big rip can be avoided in the IHDE3 model. 

In Fig. 8 we plot the evolutions of £25, (gray, left panel), Quan (green, left panel), Werfde (gray, right 
panel) and weff tor (orange, right panel) along with the scale factor a for the IHDE3 model. From the left 
panel, we see that, at the 68.3% CL O,;, — 1 is present in the future due to a negative I (see Table I), but 
the interaction can prevent €25, from approaching 1 at the large a region at the 95.4% CL. The right panel 
shows that the stable point (corresponds to Weffde = Wef f,tot = —1) can be accomplished when a > 1000 at 
the 95.4% CL. Thus, the IHDE3 model may help to avoid the big rip singularity. 

It would be worthy further investigating what kind of samples in our MC analysis may achieve the 
de Sitter solution. Of 10’ samples generated for the IHDE3 model there are about 4.7 million samples 
satisfying Ay? < 4. Within them, we find about 18,000 samples satisfying T > 0, and finally 709 samples 
having de Sitter solution in the far future. In Fig. 9, these samples are plotted in the a— (top panels) and 
I-6 (bottom panels) planes. Clearly, we see that to achieve a de Sitter solution, a large $ is required (all 
these samples have 8 > 1.5). The bottom right panel of the figure shows that, at the 68.3% (Ay? < 1) CL, 


due to a negative T, the big rip cannot be avoided. 
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FIG. 9: Ay? = 1 (olive dotted) and Ay? = 4 (green dashed) contours in the a (top panels) and I—8 (bottom panels) 
planes. Orange dots are the 709 samples which can achieve a de Sitter solution in the future and thus avoid the big rip. 
In the left panels the whole contours are plotted, while the right panels, as the zoom in of the dashed orange region of 


the left panels, show these samples in detail. 8 = 1.5 is plotted in black dotted lines. 
V. CONCLUDING REMARKS 


We have investigated the IHDE models with Q œ pes Three IHDE models, including the IHDEI 
model with Q = ['Hopgm, the IHDE2 model with Q = l'Hopaz, and the IHDE3 model with a and £ running 
freely, were investigated. 

By using the Union2.1+BAO+CMB-+AHp data, we placed cosmological constraints on these models. We 
found that a negative T, i.e., energy flow from dark matter to dark energy, is slightly favored by the data, 


although the HDE model with T = 0 still lies in the 95.4% CL region. For all the models, we get c < 1 at 
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the 95.4% CL. The values of c in the IHDE models are smaller than that in the HDE model. 

We showed that the interaction has different properties in the three models. Compared with the cosmic 
expansion, the effect of interaction on the evolution of pg, and pj, is smaller. We also put a constraint on the 
total amount of the energy change in the dark matter component. At the 95.4% CL, the increment amount 
of the dark matter is constrained to be less than 9%, while the decrement amount is constrained to be less 
than 15%. We found that the constraint basically does not depend on the forms of the interaction term. 

Furthermore, we discussed the fate of the universe by investigating the dynamical equations of the IHDE 
models. We find that the equations may give a de Sitter solution at z — —1, with the effective equation of 
state for dark energy and dark matter being —1. When confronted with data, we show that this solution 
cannot be accomplished in the IHDE1 and IHDE2 models. Rather, in these two models since the data favor 
a smaller c, the big rip is even more severe than the HDE model. In the IHDE3 model, we show that such a 
solution can be achieved for a large 6, and the big rip may be avoided at the 95.4% CL. 

Finally, there are still some issues not covered in our paper, i.e., the precise condition for obtaining the 


de Sitter solution, the coincidence problem, and so on. These issues all deserve further investigations. 
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